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ABSTRACT 

A robust detection of the tidally induced intrinsic alignments of the late-type spiral 
galaxies with high statistical significance is reported. From the spectroscopic galaxy 
sample of SDSS DR7 compiled by Huertas-Company et al. which lists each galaxy's 
probabilities of being in five Hubble types, P(E), P(E11), P(S0), P(Sab), P(Scd), we 
select the nearby large late-type spiral galaxies which have redshifts of < z < 0.02, 
probabilities of P(Scd) > 0.5 and angular sizes of D > 7.92 arcsec. The spin axes 
of the selected nearby large late-type spiral galaxies are determined up to the two- 
fold ambiguity with the help of the circular thin-disk approximation and their spatial 
correlations are measured as a function of the separation distance r. A clear signal of the 
intrinsic correlation as high as 3.4a and 2. 4a is found at the separation distance of r ~ 
1 fa Mpc and r ~ 2 /i -1 Mpc, respectively. The comparison of this observational results 
with the analytic model based on the tidal torque theory reveals that the spin correlation 
function for the late-type spiral galaxies follow the quadratic scaling of the linear density 
correlation and that the intrinsic correlations of the galaxy spin axes are stronger than 
that of the underlying dark halos. We investigate a local density dependence of the 
galaxy spin correlations and found that the correlations are stronger for the galaxies 
located in dense regions having more than 10 neighbors within 2fo~ 1 Mpc. We also 
attempt to measure a luminosity dependence of the galaxy spin correlations, but find 
that it is impossible with our magnitude-split samples to disentangle a luminosity from 
a redshift dependence. We provide the physical explanations for these observational 
results and also discuss the effects of possible residual systematics on the results. 

Subject headings: cosmology:theory — methods:statistical — large-scale structure of 
universe 



1. INTRODUCTION 



The building blocks of the large-scale structure in the universe are the galactic dark halos that 
are believed to underlie the luminous galaxies. Although the galactic dark halos are not directly 
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observable, their intrinsic properties can be inferred and estimated from the galaxy observables 
such as luminosity, color, morphology and spatial correlations. For instance, how luminous a given 
galaxy is related to how massive its underlying dark halo is. The colors and morphologies of the 
observed galaxies reflect the formation epochs and assembly history of their galactic halos. Finding 
the hidden relations between the galaxy observables and intrinsic properties of the underlying dark 
halos is one of the most important tasks in the field of the large-scale structure. 

The directions of the galaxy angular momentum (i.e., the galaxy spin axes) are one of those 
galaxy observables that are believed to have a direct connection wit h the intrinsic properties of 
the underlying dark halos. In the standard linear tidal tor que theory ([Peebles! 1 19691 : iDoroshkevich 
197Q : IWhitdll984l ; lDubinskilll992l ; ICatelan Theunslll996l ). the tidal torques from the surr ounding 
matt er originate the angular momentum of the proto-galactic halos in the linear regime (see ISchafer 
2009, for a recent review). In consequence, the directions of the angular momentum of the proto- 
galactic halos come to be al i gned with the interm ediate principal axes of the local tidal tensors 



(Lee k, Pen 



2000, 



2001 



2002 



Porciani et al.ll2002l ). As the tidal fields are spatially correlated, the 



correlations between the spin axes of the proto-galactic halos and the principal axes of the local 
tidal tensors would induce the spatial correlations of the spin axes of the proto-galactic halos. A 
detailed analysis have shown that the spatial correlations of the spin axes of the galactic halos would 
follow a quadratic scaling with the two poi nt correlation function of the linear den s ity field provided 



that the linear tidal to rque theory is valid (|Pen et al 



Crittenden et al 



2001 



2000 



Sugerman et al 



2000 



Lee Pen 



2001 



)■ 



A critical question is whether or not the tidally induced initial correlations between the spin 
axes of the galactic halos have been retained to a detectable level and how strong the retained 
correlations are. If the intrinsic correlations of the spiral galaxies are detected and if we can model 
them well with the tidal torque theory, it would make it plau s ible t o reconstruct the initial tidal 
fields from the observable spin fields as proposed bv lLee &: Pen! ([2000). Furthermore, it would have 
a direct impact on the weak lensing communit y where the intrinsic alignments of the spiral (or 



blue) galaxies are often assumed to be zero (e.g. iJoachimi et alJl2010l ; iMandelbaum et al. 



20101 ) 



Plenty of observational efforts have so far been made to address this issue. iPen et al.l (J2000) 
measured the spatial correlations of the spin axes of the nearby spiral galaxies from the Tully 
catalog and reported a tentative detection of 2a signals. They claimed that the observed trend of 
the galaxy s pin correlati o ns is consistent with the quadratic scaling of the linear density two-point 
correlations. iLee Sz Pen! (|2002l ) searched directly for the correlations between the spin axes of the 
nearby spiral galaxies and the in termediate principal axes of the tidal fields using the data from 
the Point Source Galaxy catalog (jSaunders et al.ll2000l ). 



Brown et al.l (|2002l ) demons trated that the intrin sic shape correlations of the elliptical galaxies 



from the SuperCOSMOS survey (jHamblv et al.ll200ll ) are consistent with the quadratic scaling with 



the linear density correlations under the assumption that the projected minor axes of the observed 



elliptical galaxies are orthogonal to the projected spin axes of their underlying halos. iNavarro et. al 
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(|2004l ) noted that the distribution of the inclination angles of the nearby disk galaxies relative to the 



Surjergalactic plane is consistent with the prediction of the linear tidal torque theory. iTruiillo et. al 



(|2009 j claimed that the observed alignments of the g alaxies in the vic i nity o f large voids are well 
explained by the tidally induced intrinsic correlations. iLee fc Erdogdul (|2007l ) detected a 2a signal 
of the correlations between the spin axes of the nearby spiral galaxies from the Tully catalog and 
the intermediate prin cipal axes of the local tidal field constructed from the 2MASS redshift survey 
(jErdogdu et al J 12009 . and references therein). 



Lee fc Penl (|2007l ) measured the two dimensional projected spin correlations of the blue galaxies 



from the S loan survey (lYork et al.ll2000l ) and found a 3a signal at one single distance bin of r ~ 
1 h Mpc. iPaz et al.l ( 20081 ) found the correlations between the galaxy angular momentum vectors 
and the large-scale structures from the Sloan survey and claim ed that their r esults are in qualitative 
agreement with the predictions of the tidal torque theory. ISlosar et al.1 (12009) rep orted a first 



detection of the chiral correlations of nearby galaxies. Recently, Uones et al.l (|2010i ) found that 
the spin axes of the spiral galaxies in the filaments tend to be aligned with the major axes of the 
filaments and concluded that the observed spin alignments of the filament galaxies provide a fossil 
evidence for the tidally-induced intrinsic correlations. 

In spite of these observational evidences, it is still inconclusive whether or not the present 
galaxies still retain its initially memory of the tidally induced intrinsic spin correlations for the 
following two reasons. The first reason is that the reported correlation signals are not strong 



en ough to be confirmed as true s ignals . Fo r instance, the g a 



bv lPen et al.l (|200(j ) . ILee fc Penl (|200lh and lLee fc Erdoedu 



levels. The samples on which the analyses of lNavarro et. al. 



axy spin correlation signals reported 
(2007) are significant only at 2 — 2.5a 
( 20041 ) and lTruiillo et. al.l feood ) were 



based are so small that their conclusions suffer from poor number statistics. Furthermore, there 
have rep orted some counter evid ences against the existence of the intrinsic galaxy alignments. For 
example, lArval &: Saurerl (|2005l ) searched for the alignments of the spin vectors of the galaxies in 
Abell clusters but failed in finding any signals. 

As for the stro ng correlations of the projected major axes of the elliptical galaxies measured by 
Brown et al.l (|2002l ). although it is good number statistics, the signals suffer from large uncertain- 
ties associated with regarding the projected major axes of the elliptical galaxies as the directions 
orthogonal to the spin axes of the underlying dark halos. True as it is that the minor axes of 
the elliptic al galaxies should be more or less aligned w ith the spin axes of their underlying dark 



halos (e.g.. iBailin et al 



2005 



Bailin &; Steinmeta 120051 ) . their spatial correlations are not so good 



tracers of the linear tidal fields as the spin correlations of the spiral galaxies, since the shapes of 
the elliptical galaxies are apt to be affected by the sub sequent processes such as the anisotropic 
merging/infall, galaxy-galaxy interactions and etc (e.g.. Fuller et al.lll999l ). Thus, the best targets 
for the measurement of the tidally-induced intrinsic alignments are the spiral galaxies and their 
spin orientations. 



The second reason comes from the concerns about spurious signals that could be produced by 
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systematics. There are three major sources for the systematics in the measurement of the intrinsic 
correlations of the galaxy spin axes: the weak gravitational lensing effect, presence of thick central 
bulges and inaccurate measurements of the spin axes in case that the galaxies have small angular 
sizes. It has been known that the extrinsic alignments of the galaxy shapes due to the weak gravita- 



for ellipticity-ellipticitv correlations of the elliptical galaxies) at redshifts z > 0.1 (e.g., Pen et al. 


2000; 


Croft & Metzlcr 200ol: Irleavens et al. 2000: Catelan et al. 


2001: Crittenden et al. 200ll. I2OO2I: 


•Tine 


2002:lMackev et al. 


2002l:lHevmans h Heavens! hood: Hirata & Seliakhooilffirata et aDbooi 


Mandelbaum et al. 2006 


Hirata et all 2007: Joachimi & Bridle 


2010). The presence of thick cen- 



tral bulges in the spiral galaxies could cause significant systematics since it invalidates the circular 
thin-disk approximation which is almost exclusively used to measure the orientations of the galaxy 
spin axes. In case that the galaxies have small angular sizes, their position angles and axial ratios 
would be difficult to measure accurately, which would propagate into the systematic errors in the 
measurement of the galaxy spin axes. 

The errors caused by the above three systematics in the measurement of the intrinsic correla- 
tions can be minimized if one selects only the large late-type spiral galaxies (of Hubble types Scd) 
observed at lowest redshifts. Using only low-z galaxies, one can reduce the weak gravitational lens- 
ing effect to a completely negligible level. Selecting only large late-type spiral galaxies which have 
the smallest central bulges, one can guarantee that the measurements of the spin axes through the 
circular thin disk approximation are reliable. To select such galaxies from the observational data, 
however, it is necessary to have in formation on the spectroscopic redshifts and Hubble types of 
the observed galaxies. Very recently, iHuertas-Company et all (120101 ) have released a catalog of the 
galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), in which such information 
are all available. Our goal here is to measure the true intrinsic correlations of the nearby large 
late-type spiral galaxies from this catalog and to study their behaviors and their dependence on 
the local density and luminosity. 

The organization of this Paper is as follows. In §2.1, the sample of the nearby large late-type 
galaxies selected from the SDSS DR7 is described and the number distributions of the selected 
galaxies as a function of their magnitude and local density are derived. In §2.2, it is explained 
how the spin axes of the selected galaxies are determined up to two-fold degeneracy with the help 
of the circular thin disk approximation. In §3.1, the correlations between the spin axes of the 
selected nearby large Scd galaxies are measured and the bootstrap error analysis is presented. In 
§3.2 the dependence of the galaxy spin correlations on local density and luminosity is investigated. 
In §3.3 the possible residual systematics are discussed and its effect is examined. In §4, the obser- 
vational results are compared with the analytic models based on the tidal torque theory and the 
bestfit parameters to quantify the strengths of the correlations are determined. In §5, the results 
are summarized and a conclusion is drawn. 
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PROPERTIES OF THE LATE-TYPE SPIRAL GALAXIES FROM SDSS DR7 



2.1. Selection of the Nearby Large Scd Galaxies 



The galaxy catalog compiled by iHuertas-Companv et al. (2010) consists of a to tal of 698420 
galaxies in redshift range of < z < 0.16 from the SDSS DR7 (jAbazaiian et al.ll2009l ). It provides 
information on spectroscopic redshift (z), right ascension (a), declination (<5) and the probabilities 
of being in the five morphological classes (E, Ell, SO, Sab, Scd) for each galaxy. Through private 
communication with M. Huestra-Company, we obtain additional information on the major and 
minor axes, position angle (i?p) and r-band model magnitude (m r ) of each galaxy in the catalog. 



To minimize the systematics in the measurement of the spatial correlations of the galaxy spin 
axes, we select only those galaxies which have redshifts in the range of < z < 0.02, probability 
of being in the morphological class of Scd higher than 0.5, and angular size (given as the major 
axis length) D larger than 7.92 arcsec (corresponding to 20 pixels in the SDSS frames). A total 
of 4065 galaxies from the catalog are found to satisfy these three conditions. Figure CD plots the 
number distribution of the galaxies from the catalog as a function of its probability of being in the 
Scd morphological class, P(Scd). The results shown in each panel is drawn under the condition 
of < z < 0.02 (top left panel); < z < 0.02 and P(E) < P(S0) < P(Sab) < P(Scd) (top 
right panel); < z < 0.02 and P(Scd) > 0.5 (bottom left panel); < z < 0.02, P(Scd) > 0.5 and 
D > 7.92 arcsec (bottom right panel). Using information on z, a and 5, we d etermine the comovin g 
distance to each galaxy in unit of h M.pc assuming a WMAP7 cosmology ( Komatsu et al.ll2010l ). 



To determine the local density of the selected galaxies, we first construct a volume- limited 
sample of the galaxies (regardless of their types) in the same redshift range from the SDSS spectro- 
scopic data. Basically, the volume-limited sample includes only those SDSS galaxies whose apparent 
r-band magnitudes would exceed a given flux limit, m r>c , if placed at z = 0.02. Figure [3] plots the 
fraction of the SDSS galaxies included in a volume-limited sample as a function of m rc . When the 
SDSS flux limit value of m rjC = 15.2 is applied, a total of 12273 galaxies are found to belong to 
the constru cted volume-limited samp l e. No w, for each selected nearby Scd galaxy from the catalog 
provided bv lHuertas-Companv et al.l (|2010l ). we count the number of its neighbor galaxies from the 
volume-limited sample whose separation distance from the given selected Scd galaxy is less than 
r s = 2/i _1 Mpc. Figure [2] plots the number distribution of the selected nearby large Scd galaxies 
as a function of the neighbor galaxies located within separation of r s . The distribution has its 
maximum at N ng ~ 10 and only a small fraction of the selected Scd galaxies have indeed more than 
100 neighbor galaxies within r s . It indicates that the majority of the selected nearby Scd galaxies 
are field galaxies, not belonging to galaxy clusters. 

With the measured comoving distance to each galaxy, we convert each galaxy's apparent 
magnitude m r to th e absolute magnitud e M r by using the system of the inverse hyperbolic sine 
magnitudes (asinh) (jLupton et al.1 Il999i ) . Figure |4] plots the number distribution of the selected 
nearby large Scd galaxies as a function of M r . Table [T] lists the redshift range, the total number, 
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mean absolute r-band magnitude and mean number of the neighbors within 2h Mpc, averaged 
over all the selected nearby large Scd galaxies. 



2.2. Determination of the Galaxy Spin Axes 



To determine the spin axes of each selected galaxy, we use information from the SDSS imaging 
data on galaxy's axial ratio q and position angle $p which were extracted by the SDSS team 
from the coefficients of the Fourier expansion of the 25 magnitudes per square arcsecond isophote 
measured in the r- band. Placeholders for the errors on each of these quantities were not available 



from the website (http://www.sdss.or g/dr7). We calculate the inclination angle, £, of each selected 



galaxy as (jHavnes Sz Giovaneililll984l ) 



cos 2 £ 



1-p 2 



(1) 



where p is the intrinsic flatness parameter that depends on th e galaxy morphological typ e. Here, 
we adopt the value of p = 0.10 for the Scd galaxies, given bv lHavnes Sz Giovanellil ()1984l ). If the 
inclination angle of a given galaxy in the selected sample is less than the intrinsic flatness parameter 
(i.e., q < p), then we set it at tt/2 (see §3.3 for discussion on the systematics related to the intrinsic 
flatness parameter). 

With the help of the circular thin-disk approximation and given information on the inclination 
angle £ and position angle #p, we determine the unit spin vector, L, of each selected nearby Scd 
gal axy in a local spheric al-polar coordinate system up to sign-ambiguity of the radial component 
as (|Lee fc ErdogdJboO^ 



L r 

U 

La, 



± cos £, 

(l-cos 2 1/2 sintfp, 
(1 



cos 2 £) 1//2 cos 'dp 



(2) 
(3) 
(4) 



where L r , Lg, L^ correspond to the radial, polar and azimuthal component of L, respectively. This 
sign-ambiguity in the radial components of L is due to the fact that one cannot determi ne whether 



the r otation of a given galaxy upon its symmetry axis is clock-wise or counter-clock wise (|Pen et al 

2oooh . 



Using information on the equatorial coordinates, (a, 5), we determine the unit spin vector of 
each selected Scd galaxy in the Cartesian coordinate system up to two-fold degeneracy as 



L a \ = L r sin 6 cos <f> + Lg cos 9 cos (f> — sin i 
L a 2 = L r sin 6 sin <fi + Lg cos 6 sin 4> + L^ cos <; 
L a 3 = L r cos 9 — Lg sin 6, 



(5) 
(6) 
(7) 



-7- 



Lfri = — L r sin 9 cos + Lg cos cos — Lysine 
= — -Cr sin sin + Lg cos sin <f> + cos <; 
= — L r cos6 — Lg sin 6, 



(8) 
(9) 
(10) 



where 9 = ir/2 — 5 and <f> = a. Finally, to each selected nearby Scd galaxy, we assign a set o f two 
unit spin vectors, L a and L&, which differ from each other by the sign of L r (jPen et al.ll2000l ). 



3. INTRINSIC SPIN CORRELATIONS OF THE LATE-TYPE SPIRAL 

GALAXIES 

3.1. Measurement of the Galaxy Spin Correlations 



The spatial correlation of the galaxy spin axes is defined bv lPen et al.l (|2000l ) as 



r?(r) = (jL(x).L(x + r)| 2 )-i (11) 

Here the ensemble average is taken over those galaxy pairs whose separation distance is in a 
differential range of [r, r + dr] . The value of 1/3 is subtracted since it is the value of the ensemble 
average when there is no correlation. In practice, however, due t o the two-fold de generacy in the 



determination of L, the correlation r/(r) can be measured only as (jPen et al.ll2000l ) 



V (r) = \ (<|L a • t'f) + (|L a • L' 6 | 2 ) + <|L 6 ■ t'f) + (|L 6 • t' b \ 2 )) - \ (12) 

where L' represents the unit spin vector measured at x + r. 

For all pairs of the galaxies, we calculate the squares of the dot products of their unit spin 
vectors, taking into account the two- fold degeneracy. Binning the separation distances r, we take 
the ensemble average over those pairs whose separation distances belong to a certain distance bin 
and subtract 1/3 to determine rj(r) in accordance with Equation (|12|) . Although we try to reduce 
systematic errors as much as possible through constraining the redshift-range, morphology and 
angular size, there could be some residual systematics. To sort out possible residual systematics, 
we perform a bootstrap error analysis. The spin correlation functions, rj(r), are remeasured from 
each of the 1000 bootstrap resamples which is constructed through shuffling randomly the positions 
of the selected galaxies. The bootstrap errors, erf,, are calculated as = {{rj — fjf,) 2 ) 1 / 2 where the 
ensemble average is taken over the 1000 resamples and fjb represents the bootstrap mean. If there 
were no residual systematics in the measurement of the galaxy spin axes, then the bootstrap mean 
value would be very close to zero. The degree of the deviation of the bootstrap mean from zero 
would indicate the level of the residual systematics. 

Figure [5] plots the observed galaxy spin correlations for four different cases of the galaxy 
angular size cut (D c = 0.00, 1.98, 3.96 and 7.92 arcsec in the top-left, top-right, bottom-left and 
bottom-right panel, respectively). In each panel the error bars represent one standard deviation, 



- 8 - 



a^, in the measurement of rj(r). The thin solid and dashed line in Figure [5] represents the bootstrap 
mean fjf, and bootstrap errors, 07,, respectively, while the thin dotted line corresponds to the zero 
signal. As it can be seen, a clear signal as strong as 3.4<7b is detected for all four cases of D c at the 
second radial bin corresponding to the distance of r ~ 1.25/i Mpc. For the cases of D c = and 
1.98 arcsec, the third radial bins corresponding to r ~ 2/i~ 1 Mpc also exhibit a 3<7b signal while for 
the other two cases of D c = 3.96 and 7.92, the correlation at r = 2/i~ 1 Mpc diminishes down to 
2.4CT;, level. 

As mentioned in §1, for the galaxies having small angular sizes the spin axes are difficult 
to measure accurately. As it can be seen from Figure [5j for the two cases of D c = and 1.98 
arcsec where the small angular size galaxies are included, the mean bootstrap values (thin solid 
line) show non-negligible deviations from zero (larger than 10 -4 ). Thus, the 3<Jb signals detected 
at r ~ 2/i~ 1 Mpc for these two cases are likely to be contaminated by the systematics. Figured] 
shows more clearly how the averaged bootstrap mean value changes with D c . As it can be seen, 
the bootstrap mean value decreases as D c increases and drops below 10 -4 (negligible level) when 
D c reaches up to 7.92 arcsec, which justify our choice of D c = 7.92 arcsec for this analysis. 

Since the correlation function, r/(r), is computed as a function of the three dimensional sep- 
aration, there may be some cross-correlations at different radii. If the cross-correlations between 
the second and third radial bins are significant, then the signal detection for the second radial 
bin would have lower significance. To quantify the amount of cross-correlations, we compute the 
covariance matrix, Cy, that is defined 

Cij = ((Vi ~ Vbi)(rjj ~ Vbj)) (13) 

where the ensemble average is taken over 1000 bootstrap resamples, r\i is the correlation at the i-th 
bin from a bootstrap resample, and %j is the bootstrap mean at the i-th bin. If the cross-correlations 
between the second and third radial bins are significant, then the off-diagonal component, C23, 
would be comparable in magnitude to the diagonal component, Cii- It is found that C22 = 
3.18 x 10" 6 while C 2 3 = -0.17 x 10~ 7 . This implies that the cross-correlations between the second 
and third radial bins are order of magnitude smaller (and even negative) than the correlations at 
the second radial bins, which proves that the signal detected at the second radial bin is robust. 



3.2. Dependence on Local Number Density and Luminosity 

Now that we have detected non-zero galaxy spin correlations from the SDSS data, we would like 
to investigate whether or not r](r) depends on the local number density N ng of neighbor galaixses. 
Applying the number density cut, N ngiC = 10, to the sample, we divide the selected Scd galaxies 
into two subsamples each of which consists of those nearby large Scd galaxies with N ng < N ngfi 
and N ng > N ngtC , respectively. Then, we measure r/(r) from each subsample separately, the results 
of which are plotted in the top and bottom panel of Figure [71 respectively. 
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As it can be seen, for the subsample with N ng > 10, we find a correlation signal as significant 
as 3.1<7{, signal at r ~ 1.25 /i _1 Mpc, while for the subsample with N ng < 10, no correlation signal is 
detected. The bootstrap mean values are quite close to zero for both of the cases, which indicates 
that the detected signal is not spurious. This result shows that those Scd galaxies located in denser 
regions tend to have stronger spin correlations. A possible physical explanation is that the galaxies 
located in denser regions experience so stronger tidal forces from the surrounding neighbor mass 
distribution that the tidally-induced intrinsic correlations of their sp in axes are better retained, 
which is consistent with the previous results of lLee fc Erdogdul (|2007l ). 



Using the galaxy pairs in which one galaxy has N ng > N ng ^ c and the other has N ng < N ng ^ c , 
we also measure the cross-correlations of the spin axes between the two subsamples, the result of 
which is plotted in Figure As one can see, the result is consistent with zero cross-correlation, 
which reveals that the spin axes of the galaxies in less dense regions are not correlated with those 
of the galaxies in denser regions. 

To investigate how the intrinsic correlations depend on the luminosity, we take the median 
r-band absolute magnitude, —16.67, as the threshold M r ^ c , and divide the selected galaxies into 
two subsamples each of which consists of those nearby large Scd galaxies with M r > —16.67 and 
M r < —16.67, respectively. Then, we measure rj(r) from each subsample separately, which are 
plotted in the top and bottom panel of Figure EB respectively. As it can be seen, the correlation 
function, 7](r), exhibits a 4.9<76 peak at r ~ 1.25/i _1 Mpc and a 2.4<7& peak at r ~ 2/i -1 Mpc in the 
fainter and brighter galaxy sample, respectively. 

This result, however, cannot be interpreted as a clear luminosity dependence of the galaxy 
spin correlations. In fact, the peaks of ry(r) detected at 1.25/i _1 Mpc and 2/i _1 Mpc (for the fainter 
and brighter sample, respectively) turn out to correspond to a similar angular separation of galaxy 
pairs. It indicates the existence of significant residual systematics, which is also manifest from the 
relatively high degree of the deviation of the bootstrap mean from zero for the fainter sample as 
shown in the top panel of Figure [9j Furthermore, the redshift distributions of the two subsamples 
are found to be largely different: the brighter (fainter) sample is biased to relatively higher (lower) 
redshifts, as shown in Figure [TUl The mean and median redshifts of the two subsamples are also 
listed in Table [2l Given that there should exist significant systematics and that it is impossible to 
disentangle a luminosity from a redshift dependence of the galaxy spin correlations, we admit that a 
luminosity dependence of the galaxy spin correlations cannot be measured with our magnitude-split 
samples. 

In §3.3, we discuss fully the possible sources of this residual systematic errors involved in the 
measurement of the spin correlations of the fainter galaxies. Before moving on to the discussion on 
systematics, we also measure the spin cross correlations between the fainter and brighter samples, 
the result of which is plotted in Figure [TTJ As can be seen, a strong anti-correlation of the spin 
axes as significant as 3<7b is found at r ~ 10/i _1 Mpc. In other words, the spin axes of the fainter 
galaxies tend to be orthogonal to that of the brighter galaxies separated by the distance of ~ 
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10 h 1 Mpc. Given the numerical and observational results from th e previous works ( Bailin et al 



20051 ; iBailin &: Steinmeta 120051 ; iJones et al.l [201CJ; lHahn et al.l l201Cll ) , we explain this observational 



result as follows: The spin axes of the fainter galaxies are aligned with the intermediate principal 
axes of the local tidal fields as predicted by the linear tidal torque theory. Whereas for the brighter 
galaxies which are usually located in cosmic filaments, their spin axes tend to be alig ned with the 



longe st axes of the local filaments, i.e., the minor principal axes of the local tidal fields ( Jones et al 



2010), as revealed by the recent numerical results from hydrodynamic simulations (jHahn et al 



2010 ). In consequence, the two subsamples have strong anti spin-correlations between each other. 



The separation distance of 10 h 1 Mpc at which the strong anti-correlations are detected should be 
the typical separation between the faint isolated and bright wall galaxies. 



3.3. Discussion on Systematics 

For the accurate measurement of the galaxy spin correlations, we have tried to eliminate 
systematic errors as much as we can by constraining the redshift, morphological type and angular 
size of the SDSS galaxies. The small bootstrap errors shown in the bottom right panel of Figure 
[5] suggests that no severe systematics be existent in our constrained sample, and thus that our 
detection of the spin correlations be robust. Nevertheless, it has to be admitted that our results 
may not be completely free from systematics. 

A possible source of residual systematics could be related to the constant value of the intrinsic 



flatne ss parameter, p = 0.1 in Equation ([TJ that we adopt in accordance with lHavnes Sz Giovanelli 



(| 19841 ) to accou nt for the finite thickness o f the Scd galaxies. Given that this value of p = 0.1 



was obtained by lHaynes Sz Giovanellil (119841 ) from a sample of HI 21cm observations and that the 
shape of the neutral hydrogen distribution does not necessarily coincide with the optical shape of a 
galaxy, one may suspect that a different value of p should be used for the galaxies from the optical 
SDSS survey. An ideal way would be to look at the individual images of the optical shapes of the 
SDSS galaxies and to search for the optimal value to p for the selected Scd galaxies. This task, 
however, would be extremely time consuming and thus should be beyond the scope of this paper. 

Instead, we inspect h o w the final results would change as we use different values to p in Equation 



([I]). lHaynes Giovanellil (j 19841 ) explained that unless a galaxy has a well developed bulge like the 
early type spirals, the intrinsic flatness parameter must have the value in range of 0.1 < p < 0.15. 
Since we select only Scd galaxies which optical shapes must have the thinnest bulges, we expect 
that the value of p would not differ significantly from 0.1. Figure [T2l plots the same as Figure [5] but 
for four different cases of the intrinsic flatness parameter (p = 0.0 0.05 0.1 and 0.15). As one can 
see, the final results change less than 10% as the value of p changes from 0. to 0.15, which justify 
our choice of p = 0.1. 

Another possible source of residual systematics is the limitation of the circular disk approx- 
imation. The morphology class of Scd from the SDSS catalog includes the irregular galaxies 
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(jHuertas-Companv et al.ll2010l ). for which the circular thin disk approximation would be invalid. 
But, through private communication with M.Huertas-Company, we have learned that the fraction 
of the included irregular galaxies must be small since the irregular galaxies are usually fainter while 
the targets of the spectroscopic SDSS sample are brighter. The relatively larger bootstrap mean 
values of rj(r) for the fainter samples shown in the top panel of Figure [9] are most likely to be 
caused by the inaccurate measurement of the spin axes of the irregular galaxies that are included 
in the subsample of the fainter galaxies {M r > —16.67). An ideal way should be to look at the 
individual images of the fainter galaxies and sort out the irregular galaxies from the fainter sample. 
However, it would be also extremely time consuming and thus should be beyond the scope of this 
paper. Here, we just state explicitly that the residual systematics could cause spurious signals of 
correlations for the fainter Scd galaxy sample, which is most likely to be due to the limitation of 
the circular thin disk approximation and that a better algorithm for the measurement of the spin 
axes of the irregular galaxies would be desirable. 



4. COMPARISON WITH THE THEORETICAL MODEL 



4.1. A Brief Review of the Analytic Model 

According to the linear tidal torque theory, the angular momentum vector of a proto-galactic 



the local tidal shear tensor T = (T«) 


Doroshkevich 


1970; 
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1984; 


Catelan & Theunsl 


1996; 


Lee k PenlboOO. 
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Crittenden et al. 


2001 


Lee & Pen 
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Porciani et al. 


2002; 


Lee & Erdoedu 


2007: Lee & Pen 


2008) 



Li oc ei jk TjiIi k , (14) 

where is the fully anti-symmetric tensor. Equations (|14p implies that the spin axes of the 
proto-galaxies (i.e., the direction of the proto-galaxy angular momentum vect or) are not random 
but correlated with the intermediate principal axes of the local tidal field tensors (jCatelan &: Theuns 
19961 : iLee &: Penll2000l ). Assuming that the initially generated c orrelations betwee n L and T have 



been retained to some non-negligible degree till present epoch, Lee Sz Pen (j200ol ) have proposed 
the following formula to quantify the strength of the correlations: 



(LiLj) 



1 + a 



'j > 



(15) 



where L is the unit spin vector of a given galaxy, T is the unit traceless tidal shear tensor smoothed 
on the galactic scale, and a is a free parameter to quantify the strength of the correlations between 
L and T. 



A direct application of Equation (|15p to the spatial correlations between the sp in axes of 
neighbor galaxies have yielded the following formula for the galaxy spin-spin correlations (jPen et al. 
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2000; iLee fc Penll200ll : ICrittenden et aUEOOll ) 



r/(r) 



a 2 £ 2 (r;R) 

6 e(0;R)' 



(16) 



where £(r) is the two-po int correlations of th e linear density field smoothed on the galactic scale 
(R « 1 /i _1 Mpc) (e.g..see lBardeen et al.lll98a ). According to Equation ([16]) , the spatial correlations 
between the galaxy spin axes would follow a quadratic scaling of the density correlation, decreasing 
rapidly to zero as r increases. We fit the observations results obtained in §3.1 and 3.2 to Equation 
(|16p to determine the best-fit correlation parameter, a and examine whether or not the observed 
spin correlations follow the quadratic scaling with the linear density correlation, £(r). 



4.2. Observed Scaling of the Galaxy Spin-Spin Correlations 



We minimize the following generalized \ 2 to find the bestfit-value of a (jHartlap et al.l 120071 ). 

l^6j n — 1) 

where C^ 1 is the inverse of Cy that is defined in Equation (|13p . r/(rj) represents the theoretical 
value calculated at the i-th distance bin, r)i is the observed correlation at the same distance bin 
(with the bootstrap means subtracted) and Nun is the number of radial bins. Note that since 
there is only one parameter, the degree of freedom for x 2 ls given as Ny, n — 1. Here, the covariance 
matrix is computed using the bootstrap resamples as defined in §2.2 and the fitting is done for 
the distance range of 0.5 < r/(/» -1 Mpc) < 50 since at larger distances the numerical flukes tend 
to contaminate the fitting. If this reduced x 2 have a value close to unity, the theoretical model, 
Equation (fT6|) . would be regarded as a good fit to the observational result. 



Figure [13] plots the best-fit model as solid line and compares it with the observational results 
(square dots). The best-fit value of a is found to be 0.25 ±0.04 and the corresponding value of \ 2 is 
fou nd to be 0.83 (see T able El). Here, the error associated with the measurement of a is determined 



as (jHartlap et all 120071 ) 



Now that the reduced x 2 is found to be close to unity and the best-fit value of a is found to be 
6a a higher than zero, it can be said that the observed spin correlations are indeed tidally induced, 
having quadratic scaling with the linear density correlations, as predicted by the linear tidal torque 
theory. 



We note here that both of the best-fit value of a are larger than that determined bv lLee Pen 



(|2008l ) through fitting of the numerical data from the Millennium Run simulations (see Table 3 
in Lee & Pen 2008). There are two important implications of this result. First of all, the spin 
orientations of the luminous galaxies may not be perfectly aligned with those of their dark halos, 
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which is consistent with the previous numerical results (e.g.. iBailin et al.ll2005l ; iBailin fc Steinmetz 
20051 : lHahn et al.ll2010l ). Second of all, the luminous galaxies have retained better the initial mem- 
ory of the tidally induced intrinsic spin-spin correlations than their dark halos. This phenomenon 
may be explained as follows. The luminous galaxies are located in the inner parts of the galac- 
tic halos while the dark matter particles a re stretched to the outer parts of the dark halos (e.g., 



Bailin Steinmetz 



2005 



Hahn et al.ll2010l ). Therefore, the spin orientations of the luminous galax- 



ies may be less vulnerable than those of their dark counterparts to the destructive nonlinear effects 
from the surroundings which would break the tidally induced intrinsic correlations. 

The fitting results for the subsamples with N ng > 10 and M r > —16.67 are also plotted in 
Figures [T4~l and [15"| respectively. The best-fit values of a and the minimum value of x 2 corresponding 
to these two cases are also listed in the second and third rows of Table El As one can see, for the 
case of N ng > 10, the best-fit value of a has the largest value, approximately 0.34. As mentioned 
in §3.2, it may be due to the stronger tidal forces from the surrounding. 



5. SUMMARY AND CONCLUSIONS 



Selecting those large nearby late-type spiral galaxies which have angular sizes larger than 7.92 
arcsec, low redshifts (0 < z < 0.02), and morphological type of Scd from the SDSS DR7, and we 
have measured the spatial correlations between their spin axes. The summary of our results is the 
following. 



• A clear signal of the intrinsic correlation as strong as 3Aaf, and lAa^ (where is the boot- 
strap error) is detected at r ~ 1.25/i~ 1 Mpc and r ~ 2.0/i -1 Mpc, respectively. The cross- 
correlations between the two radial bins are found to be insignificant. 

• The galaxies having more than ten neighbors within separation of 2 /i _1 Mpc exhibit higher 
intrinsic spin correlations (at 3.1(7;, level) than those having ten or less neighbors within the 
same separation distance. This local density dependence of the galaxy spin correlations may 
be due to the stronger tidal forces from the denser environments. 

• The observed spin correlation functions of the selected Scd galaxies follow the quadractic 
scaling with the linear density correlations, which is consistent with the prediction of the 
linear tidal torque theory. 



We have also attempted to measure a luminosity dependence of the galaxy spin correlations by 
splitting the galaxy samples according to their absolute r-band magnitudes (M r ), and noted that 
the correlation function has a peak at 1.25 /i~ 1 Mpc and 2 /i _1 Mpc for the fainter (M r > —16.67) and 
brighter sample (M r < —16.67), respectively. However, the two peaks from the two magnitude- 
split samples turn out to correspond to a similar angular separation, which hints at significant 
systematics. To make matters worse, the redshift distributions of the two magnitude-split samples 
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are found to be so different that it is impossible with our samples to disentangle a luminosity from 
a redshift dependence of the galaxy spin correlations. Therefore, we conclude that a luminosity 
dependence of the galaxy spin correlations cannot not be measured from our samples. 

The behaviors of the detected correlations are found to be consistent with the predictions of 
the tidal torque theory. But, the correlation strength turns out to be greater than predicted by the 
numerical simulations. This result indicates that the spin orientations of the luminous galaxies are 
not aligned with those of the underlying dark halos and that the luminous galaxies retain better the 
initial memory of the tidally induced spin alignments than their dark counterparts. It is physically 
explained as the spin orientations of the luminous galaxies may be less vulnerable than those of 
their dark counterparts to the destructive nonlinear effects from the surroundings which tend to 
break the tidally induced intrinsic correlations. It has to be understood fully in the future how and 
under what circumstances the spin axes of the late-type spiral galaxies have retained their original 
orientations. 

It is also worth mentioning the apparent inconsistency between our results and the previous 
studies which claimed that the two dimensional pro jec ted shapes of the b l ue SP SS galaxies show 



no correlation at higher redshifts, z ~ 0.1, (see e.g., iMandelbaum et al.l |2010| ). Given that we 
have detected 3. 4a correlation signal at z < 0.02 and that the intrinsic spin correlations of the blue 
galaxies would be stronger at higher redshifts, one may naturally expect that the shape correlations 
of the blue galaxies at z ~ 0.1 would be strong enough to affect the weak lensing signals, which 
is inconsistent with the previous results. We think of several reasons for this inconsistency: First, 
the two dimensional projected shapes of the blue galaxies may not be a good indicator of the 
tidally induced alignments; Second, there may be considerable uncertainties involved in the accurate 
measurements of the shape correlations of the blue galaxies at z ~ 0.1. As a final conclusion, our 
results will shed a new light on the study of the galaxy intrinsic alignments and the weak lensing 
analysis as well. 
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Fig. 1. — Number distribution of the nearby galaxies (0 < z < 0.02) from the SDSS DR7 as a 
function of P(Scd), the probability of being classified as Scd galaxies for four different cases: all 
nearby galaxies (top left panel); P(Scd) is largest (top right panel); P(Scd) > 0.5 (bottom left 
panel); with diameter larger than 20/i _1 Kpc and P(Scd) > 0.5 (bottom right panel). 
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Fig. 2. — Number distribution of the selected Scd galaxies from the SDSS DR7 as a function of the 
number of neighbor galaxies located within the separation distance of 2 /i _1 Mpc at < z < 0.02. 




Fig. 3. — Fraction of the galaxies belonging to the volume-limited SDSS sample as a function of 
the r-band magnitude limit, m r>c . The dashed line corresponds to the applied magnitude cut. 




Fig. 4. — Number distribution of the selected Scd galaxies from the SDSS DR7 as a function of the 
r-band absolute magnitude at < z < 0.02. 
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Fig. 5. — Correlations of the unit spin vectors of the selected Scd galaxies (diamond dots) as a 
function of the separation distance r at < z < 0.02 for the four different cases of the galaxy 
diameter cut, D c . In each panel, the thin solid line represents the mean correlation averaged 
over the 1000 bootstrap resamples, the dashed lines represent one standard deviation among the 
bootstrap resamples, and the dotted line indicates the zero correlation signal. 
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Fig. 6. — Bootstrap mean averaged over the 1000 bootstrap resamples measured at r ~ 10 h Mpc, 
as a function of D c . 
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Fig. 7. — Same as Figure[5]but using only those Scd galaxies which have more (less) than N ngtC = 10 
neighbor galaxies within separation r s = 2 TtMpc in the upper (lower) panel. 
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Fig. 8. — Same as Figure but between the Scd galaxies having more than 10 neighbors and the 
Scd galaxies having less than 10 neighbors. 
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Fig. 9. — Same as Figure [5] but using only those Scd galaxies which are fainter (brighter) than 
M rjC = —16.67 in the upper (lower) panel. 
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Fig. 10. — Redshift distribution of the fainter galaxy sample (dashed line) and bright galaxy sample 
(solid line). 
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Fig. 11. — Cross-correlations of the unit-spin vectors between the faint and the bright Scd galaxies 
at < z < 0.02. 
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Fig. 12. — Same as the bottom-right panel of Figure [5] but for four different values of the intrinsic 
flatness parameter, p. 



- 31 - 



0.0100 -~ 



0.0010 r 



0.0001 r 




10 



r [Mpc/h] 



Fig. 13. — Comparison of the spin-correlations of the selected Scd galaxies (square dots) with 
the theoretical model with the best-fit parameter (solid line). The data points (square dots) are 
obtained by subtracting the bootstrap average from the observed spin-correlation values. The errors 
represents la bootstrap scatter. 
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Fig. 14. — Same as Figure [13] but only for those Scd galaxies with N ng > 10. 
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Fig. 15.- 



Same as Figure [13] but only for those Scd galaxies with M r > —16.67. 
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Table 1. number of the selected Scd galaxies, mean redshift, mean r-band absolute magnitude, 

and mean number of neighbor galaxies within r s . 













z 


M r 




4067 


0.01 


-16.56 


27 
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Table 2. selection condition, number of the galaxies, median redshift, and the mean redshift 



condition 


N 9 


Zjn 


z. 


M r > -16.67 


2034 


0.007 


0.009 


M r < -16.67 


2033 


0.015 


0.014 
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Table 3. selection condition, number of the galaxies, best-fit value of a, reduced x 2 . 



condition 


N 9 


a 






4067 


0.25 ± 0.04 


0.83 


N ng > 10 


1945 


0.34 ±0.07 


1.10 


M r > -16.67 


2034 


0.27 ±0.04 


1.87 



